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Abstract 

The Klein-Gordon-Einstein equations of classical real scalar fields have time- 
dependent solutions (periodic in time). We show that quantum real scalar 
fields can form non- oscillating (static) solitonic objects, which are quite simi- 
lar to the solutions describing boson stars formed with classical and quantum 
complex scalar fields (the latter will be studied in this paper) . We numerically 
analyze the difference between them concerning the mass of boson stars. On 
the other hand, we suggest an interesting test (a viable process that the boson 
star may undergo in the early universe) for the formation of boson stars. That 
is, it is questioned that after a second-order phase transition (a simple toy 
model will be considered here), what is the fate of the boson star composed 
of quantum real scalar field. 
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I. INTRODUCTION 



The presence of dark matter has been estabhshed indirectly in a wide range of scale of 
the universe, from that of individual galaxies to the entire universe itself [1] . Though direct 
measurements of the nature of the dark matter have not yielded any result, speculations 
on its composition vary from baryonic to non-baryonic matter. Particles like axions and 
neutralinos are the specific targets for direct observation since the indirect measurements 
like rotation curves of spiral galaxies and others do not depend on the presence of a particular 
type of particles. One of the most promising candidates for dark matter is the boson star, 
which was discovered theoretically over thirty years ago [2] [3]. Until now, the reality of 
boson stars has been successfully applied to various plausible physical situations [4]. 

The boson star is a self-gravitating compact sohtonic object made up of bosonic fields. 
Non-interacting complex scalar fields [2] [3] were originally considered for the constituents 
composing boson stars. In this case, the resultant configurations are typically 'mini'-boson 
stars, which have small size and mass. This result originates from the following specific 
feature of boson stars; the boson star is protected from gravitational collapse by the Heisen- 
berg uncertainty principle, instead of the Pauli exclusion principle that applies to fermionic 
stars, and the characteristic length scale of the former is much smaller than that of the 
latter. It has been shown that this situation can be dramatically changed by introducing 
self-interacting complex scalar fields. The self-interaction effectively generates a repulsive 
force and the maximum mass of stable boson stars can be enhanced up to a size of the order 
of ordinary fermionic stars [5]^. 

On the other hand, there is another type of gravitationally bound solitonic object known 
as oscillating soliton star composed of classical real scalar fields [7] . In this case, the space- 
time geometry and real scalar field satisfying the Klein- Gordon- Einstein (KGE) equations 



^Analyltic evaluation of the maximum mass and higher order self-interaction effect to boson star 
configuration can be found in [6]. 
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are time-dependent (periodic in time). Such objects are of interest to study dark matter, 
since the most promising candidates for dark matter are described by real scalar fields such 
as the axion [8]. However, even though these solutions are stable with respect to some 
simple perturbations, it is still unclear whether the stability of the oscillating star can be 
maintained with general perturbations. Moreover, due to its oscillating feature, real scalar 
fields may not form primordial solitonic objects in the early universe [10] (cf. [9]), which 
are expected to play important roles in galaxy formation, the microwave background, and 
formation of protostars. 

In this paper, we report some very interesting results for the self-gravitating solitonic ob- 
ject formed with real scalar fields: there exist gravitationally bound non- oscillating (static) 
objects composed of quantum real scalar fields, instead of classical real ones. Our (zero- 
node) solitonic solution is quite similar to the solutions describing boson stars formed with 
classical [5] and quantum complex scalar fields (the latter will be studied in this paper). 
The only difference between them at the equation of motion level is just the effective coeffi- 
cient of the self-interaction term (within the validity considered in this paper, i.e. zero-node 
solutions, the semi-classical and Hartree (mean- field) approximations). Therefore, non-self- 
interacting quantum/classical complex and quantum real scalar fields coupled to gravity can 
form identical static mini-boson stars. 

On the other hand, as many viable cosmological scenarios indicate, quantum fields 
present in early stage of the universe may experience phase transitions by temperature 
changes. Thus, if the boson star forms in the early universe, it would undergo phase transi- 
tions. It might be an interesting issue to consider the consequence of such a phase transition 
for the boson star.^ What is the fate of the scalar field and what are the final products 
after the phase transition, which is presumably second order? In this paper, as a first step. 



^In [11], a boson star in an evolving cosmological background with time varying gravitational 
constant (i.e., in the context of scalar-tensor theory of gravity) has been considered. 
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we consider a simple toy model; without speculating on the detailed procedure of the phase 
transition, we introduce a plausible form of the Lagrangian for real scalar fields coupled to 
gravity, which would give the physics after phase transition, and examine that there may 
still exist boson stars. If they continue to exist, what are the differences between the boson 
stars before and after phase transition. 

In the next section, we start by studying boson stars composed of quantum complex scalar 
fields, and compare them with the boson stars formed with the classical complex scalar fields 
[5]. The boson stars composed of quantum real scalar fields are studied in Sect. 3. In Sect. 4, 
we consider boson stars composed of quantum real scalar fields that undergo a second-order 
phase transition. Disscussions and some remarks are given in Sect. 5. 



where we assume that the coupling to gauge fields is negligible. Then, the field equations of 
the scalar field are given by 



and its complex conjugate. For convenience, we introduce two real scalar fields, defined by 



II. COMPLEX QUANTUM SCALAR FIELDS 



Consider a self-interacting complex scalar field Lagrangian given by 




(1) 




(2) 




(3) 



and rewrite the field equations as 




(4) 




General expressions for solutions to Eqs. (4) can be written as 
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$1 = E (an/„.C(r)n/lle— + alU{r)li{Ylr e^^'A , (5) 

nlm 



nlm 

As quantum scalar fields, <l>i and $2 become operators 



^1 = E (anU{rUYie-^-n.^ + al^rUYiYe^^A , (6) 

nlm 

where ajj^m) ^Lm ^'^d Sni^, are the creation and annihilation operators, respectively, 
and they satisfy the commutation relations [Sj^,aa'] = [ba,ba'] = Saa'- In the following, only 
the ground state of the system, n — 1,1 — m — 0, will be considered and the condition 
ojqi = ojqi = o; is assumed. In addition, we assume that in the ground state, which is 
denoted by >= |A^$^ > (8)|A^$2 tbe number of particles created by is the same 
as that associated with b^, i.e. A^$^ = N^^ — ^- 1^ general, C('^) and r){r) are complex 
functions, but here we restrict them to real ones. To linearize the field equations (4), we 
adopt the Hartree approximation such that ($1 + $1 ~< N\{^1 + ^l)\N > $1. 

On the other hand, we also use the semi-classical approximation in Einstein's equations 

= SttG < N\fi;\N >, (7) 

where is the Einstein tensor and is the energy-momentum operator that is obtained 
by quantizing the energy-momentum tensor 

7^ = I {d^<t>*dA + d'^<t>d,<t>*) - ^gi: + \cpf + ^ |0|^ j 

The background spacetime is restricted to be static and spherically symmetric. Thus, in 
Schwarzschild coordinates, it can be written as 

ds" = -B{r)dt^ + A{r)dr^ + r^dVL\, (9) 

where ciO| denotes the metric of a two-dimensional unit sphere. 



After performing some algebraic calculations and reparametrizations in Eqs. (7) and (4), 
we obtain the equations 

A' 



B' 

ABx 



(t[ + A 



- 1 (Ti - A((7i + (72)(7i 



(7o + - H C7o + A 



'^2 

^ - 1 j (72 - A((72 + al)a2 



= 0, 
= 0, 



(10) 



where x = mr, = (AttGN)^^'^ ({r), (72 = (AnGNf^'^riir), O = u/m, A = A/47rG'm2, 
and the prime denotes the derivative with respect to x. In Eqs. (10)^, we have used the 
approximation, ^ 1. In this paper, we consider only a special type of solution given by 
(Ji — a2 — (7/-\/2. Then, the equations become 



A' If 1 



A^x x^ 
B' 1 



'Q^ \ 2 5A 4 {a'Y 



ABx 



x^ 



1 - 



5A 



B 8 A 



a" + {--{ -]a' + A 

\x 2B 2A 







(11) 

(12) 
(13) 



Note that comparing (11)-(13) with the KGE equations for classical complex scalar fields 
[5], only the coefficients of self-interaction terms in (11) and (12) are changed as A/2 
5 A/8, while the field equation (13) remains without any change. Such a change results in 
modifications of some properties of the boson stars, especially the maximum mass of stable 
boson stars. This can be seen in the following numerical analysis. 



^Note that while Eqs. (4) and (8) are invariant under a rotation in the internal space given by 
$1 $icosB + $2sin0, $2 — $isin6 + $2cos0, where G is a constant, in Eq.(lO) such 
a symmetry does not exist. This is simply because a\, 02 in Eq.(lO) are not field variables but 
(redefined) radial functions of the mode solutions. So, there is no reason maintaining the internal 
symmetry in the level of Eq.(lO). 
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Substituting an ansatz of the metric function A by 



2M {x) 



1 -1 



A{x)^ 1 



X 



(14) 



into the above equations and using the boundary conditions 



M{0) = 0, (7(0) = (7c, a'{0) = 0, and B{oo) = 1, 



(15) 



we obtain the mass of boson stars as a function of Uc for A = 100 in Fig.l. In Fig.l, it is 
shown that the maximum mass of boson stars composed of quantum complex scalar fields is 
greater than that in the case of classical complex scalar fields. This is true for other values 
of A as can be seen in Fig. 2. The lines in Fig. 2 are chosen in a way that they hold for A 1; 
in the case of quantum fields (5A/8) M^^^x ~ 0.223A^^'^ M^/m, while M^^ ^ 0.22A'^/'^M^/m 
for the classical fields (A/2) [5]. It has also to be mentioned that in Fig.l, there is a crossing 
(at (7c ~ 0.132) of the two mass-curves, and when (7c > 0.132, the mass of the boson stars 
formed by the classical complex scalar field is greater than that of the boson stars formed 
by a quantum complex scalar field. 

As mentioned above, the difference between the two cases of classical and quantum fields 
is in the coefficient of the self-interaction term in the Einstein equations (11) and (12), 
i.e. A/2 for the classical field and 5 A/8 for the quantum one. Why does the difference 
appear? In fact, the KGE equations for the classical complex scalar field considered in Ref. 
[5] become the equations for a "real" field^ after ehminating the time-dependent part. This 
is because the authors of Ref. [5] have taken the solution of the complex field which has the 
form of (p{r,t) — $(r) exp(— ia;i) with real $(r). In our case, the KGE equations given in 
(10) are for two "real" fields. Even though we chose a special form of the solution given by 
(Ti = (72 = (j/^/2 to compare with the classical case studied in [5], the resulting equations 
given by (11)- (13) already include the interaction between the two "real" fields as well as 



^The quoted term of "real" field will be used for the spatial part of a mode function, not for an 
ordinary field. 
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self-interactions of each fields, and this is the reason why the coefficients of self-interaction 
term are effectively different from each other in the quantum and classical cases, even though 
we have started with the same Largrangian (1). 

Roughly speaking, since larger self-interaction energy exerts greater repulsive force in 
the formation of boson stars, the 'quantum' effect, i.e. using the quantum field rather than 
the classical one to form a boson star, enhances the maximum mass of the boson stars. 
(In fact, when the self-interaction is introduced, density of the boson star is approximately 
proportional to A"^ so that the boson stars become dilute. However, since its radius is ap- 
proximately proportional to A^/^, the maximum mass of the boson star is to be proportional 
to A^/^. This argument is available for the case that the self- interaction energy is comparable 
to the mass energy of the field, i.e. A >> 0. See Ref. [6] for a detailed analysis.) 

Of course, the validity of our argument has to be restricted to solutions that are consid- 
ered in this paper and [5]. However, it is reasonable to study the difference between the roles 
of classical and quantum fields in the formation of boson stars. Another point that has to 
be mentioned is that the mass curves in Fig.l cross each other at Uc ~ 0.132, and the boson 
star masses for the classical and quantum cases reverse in magnitude. This phenomenon 
can not be described by the above argument. To understand such a behavior, a detailed 
analytical approach would be required. 



Now, let us consider a real scalar field for the formation of boson stars. The real massive 
scalar field Lagrangian including self-interaction is written as 



III. REAL QUANTUM SCALAR FIELDS 



(16) 



And the energy-momentum tensor and field equation are obtained as 




(17) 



(18) 
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Prom the general form of solutions of the field equation, the field operator is given by 



(19) 



nlm 



where the creation and annihilation operators, dl^i^ and dnim, satisfy the commutation re- 
lations [d]^,da'] = Saa'- Again we only consider the ground state. Note that we take the 
function $(r) as a complex function. Then, defining new real functions 



e(r) = ^ ($(r) + r(r)) , x{r) ^ ^ (*(r) - $*(r)) , 



(20) 



and taking the Hartree and semi-classical approximations, we obtain the KGE equations 
given by 

+ ^ (i - i) = (¥ + 1) + ^i) + X (p? + Pi)' + i ((p'l)' + iP2f) , 
^ - (1 - i) = (¥ - 1) (p? + pl) - f (p? + plf + i ((p'l)' + ip'^f) > 

Pi + (i + i - ^) Pi + ^ [(¥ - 1) Pi - HpI + pDpi] - 0, 

P2 + (i + i - ^) P'2 + ^ [(¥ - l) P2 - A(P? + Pi)P2] = 0, (21) 

where x = mr, pi = {AttGN)^^"^ ^{r), p2 = (47rG'A^)^^^ x(r), O = uj/m, A = A/47rGm^, and 
we use the metric (9) and the approximation N ^ 1. Taking a special form of the solution 
for P1 — P2 — c/-\/2, above equations become 



-" + (i + S-^)-' + ^[(¥-i)--A-^: 



0. 



(22) 
(23) 
(24) 



Note that the above equations can be obtained from (10) by substituting cti — (t,(72 — 0. 

Comparing Eqs. (22)-(24) with (11)-(13), wc find that the only difference is the coeflicient 
of the self-interaction term. Thus, the KGE equations (22)- (24) should have solutions that 
are much like the solutions describing the boson star composed of classical/quantum complex 
scalar fields. As mentioned above, since the coefficient 3A/4 is greater than those in the 
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classical and quantum complex scalar field cases, which are A/2 and 5A/8, respectively, 
the maximum mass of boson star formed by quantum real scalar fields should be larger 
than those of other cases. Indeed, using the ansatz (14) and the boundary conditions given 
by (15), we calculate numerically the mass of boson star as a function of Uc for various 
values of A in Fig.3. In Fig.l and Fig.2, we compare the maximum mass of boson stars 
composed of quantum real scalar fields with other cases, and confirm that the order of 
magnitude of maximum masses follows that of the coefficient of self-interaction term, i.e. 
3A/4 > 5A/8 > A/2. (The solid line in Fig.2 is the relation M^^ ^ 0.225 A^^^M^/m, which 
holds for A >> 1 for the case of quantum real scalar fields.) This result can be interpreted 
in the same way as given in Sect. 2; comparing (21) with (10), the Einstein equations in (21) 
include the term representing the interaction of "real" fields, (3A/2)p^p|, which has a larger 
coefficient than in other cases, e.g. Aerial in the Einstein equations in (10). 

IV. PHASE TRANSITION 

In this section, we consider a phase transition that the boson star can undergo in early 
stages of the universe. Our model of the system after the phase transition is a simple one 
with Lagrangian for real scalar fields given by 

>Cltt = —d^^d,^ + Im'^' - (25) 

In order to simulate the condition of phase transition the mass m in the Largrangian is 
taken to be a function of time t. Then, in the process of the phase transition, we consider 
— > — after phase transition. Such a procedure is commonly used in describing a 
second-order phase transition, and compares well to the ideas of spontaneous symmetry 
breaking which is the theoretical underpinning for such a phase transition. 

In this theory, the vacuum expectation value of the field does not vanish, but < $ >vac= 
my/X. Now, shifting the field, $ — > $ -|- m/VX, and the potential energy, V{^) — > V{^) — 

/ (4A) , we obtain a new Lagrangian written as 
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CZ^u = -\d^'^df,^ - rri'^^ - - mVX^\ (26) 

Then, the energy-momentum tensor and field equation are evaluated from the Lagrangian 
(26); 

Ti; = - i< (d"^da^ + 2m2$2 + A$4 ^ 2my/\^^) , (27) 

V^V^^ - 2m2$ - - 3mv^$2 ^ g. (28) 

Comparing (27) and (28) with (17) and (18) for the case before the phase transition, we see 
that the scalar field becomes effectively more massive, m — > -\/2m, through the phase tran- 
sition, and an additional non-linear interaction term with a coupling constant proportional 
to my/X is included in the energy-momentum tensor and field equation. Within the validity 
of the Hartree and semiclassical approximation, however, the non-linear coupling term pro- 
portional to m\/X docs not contribute to the KGE equations, i.e. < iV|$^|iV >= in the 
Einstein equations and <A/^|$|A^>$ = Oin the field equation due to the semiclassical and 
Hartree approximation, respectively. Thus, effectively the consequence of phase transition 
is just to enhance the mass of scalar fields. 

Explicitly, the KGE equations are written as 

A' 

A^3 



ls + i(i-i) = (¥ + 2)-^ + ¥-* + *^. (29) 

:fij-*(i-i) = (f-2)-^-¥-^ + ^. (30) 

+ (i + S - S?) + -4 [(¥ - 2) <^ - Aa'] = 0. (31) 

More conveniently, reparametrizing the above equations by y = \/2x, fl = 0/\/2, we obtain 



_A 
A 



ABy 



^ + ^(l-i) = (f + l)-^ + f-^ + ^, (32) 
;^-*(l-i) = (|-l)-^-f-^ + ^, (33) 



(f + S - £) + - = 0, (34) 

where prime denotes the derivative with respect to y. Note that the coefficient of the 
self-interaction term in (34) is not —A but —A/2, which is different from (13) and (24). 
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The KGE equations given in (32)-(34) are obtained by just replacing A — > A/2 in all the 
equations (22)- (24). In other words, the phase transition makes the scalar field effectively 
less self-interactive. This result apparently tells us that after phase transition the boson 
star composed of real scalar fields is to have a smaller mass than in the case before phase 
transition. 

Such a loss in the boson star mass caused by the phase transition, however, does not 
appear to be a serious shortcoming for describing the dark matter. Rather, the fact that 
the boson star can still exist after a phase transition, seems to deserve our attention. 

V. DISCUSSIONS 

In this paper, we have considered self-gravitating solitonic objects made up of quantum 
complex/real scalar fields. It has been shown that quantum complex/real scalar fields may 
compose the boson stars that are similar to that formed by classical complex scalar fields. 
The difference in the theories considered here appears only in the coefficient of the self- 
interaction term in the Einstein equations, i.e. they are A/2 for classical complex fields, 
5A/8 for quantum complex ones, and 3A/4 for quantum real ones. Numerically we have 
verified, Fig.l and Fig.2, that the maximum mass of the boson stars increases with the 
magnitude of the coefficient of the self-interaction term. This result can be understood 
as follows; after eliminating the time-dependent part the KGE equations for the quantum 
complex/real scalar fields become effectively for two "real" fields as (10) and (21). Then, 
the equations include interactions between the two "real" fields as well as self-interactions of 
each fields. On the other hand, the KGE equations for the classical complex fields considered 
in Ref. [5] become effectively for one "real" field after ehminating the time-dependent part 
and contain just the self-interaction term of the "real" field. Therefore, the coefficients of 
the resulting interaction terms in the KGE equations for the quantum complex/real scalar 
fields, which are given by (11), (12), and (22), (23), are greater than that in the case of 
the classical complex scalar fields. According to the argument that more interaction energy 
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generates effectively a larger maximum mass of the boson star, the maximum masses of the 
boson stars composed of quantum complex/real scalar fields are greater than that in the 
case of classical complex scalar fields. This argument, however, has a limitation due to the 
presence of a crossing between the mass-curves in Fig.l and Fig.2. 

In fact, the existence of such solitonic solutions does not guarantee that the bosonic 
compact objects can be formed in the universe, (cf. [12]) It has been shown that there 
is a process that is able to describe the formation of the bosonic compact objects in the 
early universe [10]. The process, which is similar to the way of describing the settling of 
coUisionless star systems, starts with collapsing due to a gravitational instability analogous to 
the Jeans instability. Then, after undergoing the so called gravitational cooling mechanism 
ejecting part of the scalar field, bosonic compact objects could be formed from the primodial 
bosonic cloud. In the case of complex scalar fields, this mechanism works to form mini- 
boson stars. However, in general, the oscillatons made up of classical real scalar fields can 
be formed only in a short dynamical time scale, but are unstable in such a state. Thus, 
without introducing additional proviso such as fragmentation of the primodial bosonic cloud, 
the classical real scalar fields would be ruled out as a candidate for the dark matter. 

In this paper, we have shown that quantum real scalar fields can form the hoson star 
rather than the oscillaton. Especially, without the self-interaction A = 0, the mini-boson 
star composed of the quantum real scalar fields becomes exactly identical to that made up 
of the classical (and quantum) complex scalar fields. Thus, considering the quantum effect, 
the real scalar field can be saved, and can be the most promising candidate for dark matter. 

As a corollary, which would be an interesting issue in the study of bosonic objects, 
the quantum complex/real scalar fields in excited states may be considered. The solitonic 
solutions, if they exist, are unlikely to be static in both (complex and real) cases. However, 
it is not obvious if such quantum fields in excited states can be ruled out completely as 
candidates for the dark matter. 

Another interesting test for the boson star composed of the quantum real scalar fields is 
the second-order phase transition. In this analysis, we have shown that the phase transition 
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makes the fields effectively more massive, m — > y/2'm, or equivalently less self-interactive, 
A — >■ A/2, and more importantly the boson stars can exist even after the phase transition. 
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FIGURES 
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FIG. 1. Boson star masses as a function of cTc when A = 100. 3A/4, 5A/8, A/2 denote the cases 
of real and complex quantum scalar fields, and classical complex ones, respectively. 
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FIG. 3. Mass of boson stars composed of real quantum scalar fields as a function of (Jc when 
A = 0,10,50,100,250. 
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